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Abstract. Let A be a (possibly unbounded) self-adjoint operator on a separable 
Hilbert space ff. Assume that a is an isolated component of the spectrum of A, 
that is, dist(o\E) = d > where E = spec(A) \ <J. Suppose that V is a bounded 
self-adjoint operator on f) such that ||V|| < d/2 and let L = A + V, Dom(L) = 
Dom(A). Denote by P the spectral projection of A associated with the spectral 
set cr and let Q be the spectral projection of L corresponding to the closed ||V||- 
neighborhood of cr. Introducing the sequence 

x n =U\- S"!~?n . ne{0}uN, 
we prove that the following bound holds: 

arcsin(||P-e||) < M„ ( 

where the estimating function M it (x), x 6 [0, , is given by 

, , 1 , , / 4k \ 1 / K ( x - H nAx)) 

M+yx) = - n t (x) arcsin — = I H — arcsin 



2 * w ' W + 2 ^ l-2x n#w ) )' 

with n t (x) = max {n I n e {0} UN,z n < x} . The bound obtained is essentially 
stronger than the previously known estimates for \\P — Q\\. Furthermore, this 
bound ensures that Q\\ < 1 and, thus, that the spectral subspaces Ran(P) 
and Ran (Q) are in the acute-angle case whenever || V|| < c+d, where 

^-2^ 4 + 32^-32 
= 16 —~ — —. = 0.454169 .... 

(7T 2 +4) 4 

Our proof of the above results is based on using the triangle inequality for the 
maximal angle between subspaces and on employing the a priori generic sin 26 
estimate for the variation of a spectral subspace. As an example, the boundedly 
perturbed quantum harmonic oscillator is discussed. 
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1. Introduction 

One of fundamental problems of operator perturbation theory is to study variation 
of the spectral subspace corresponding to a subset of the spectrum of a closed linear 
operator that is subject to a perturbation. This is an especially important issue in 
perturbation theory of self-adjoint operators. 

Assume that A is a self-adjoint operator on a separable Hilbert space Sj. It is 
well known (see, e.g., Ifl3ll ) that if V is a bounded self-adjoint perturbation of A then 
the spectrum, spec(L), of the perturbed operator L =A+V, Dom(L) = Dom(A), 
lies in the closed 1 1 V | | -neighborhood 1| v|| ( s P ec (^)) °f me spectrum of A. Hence, 
if the spectrum of A has an isolated component a separated from its complement 
E = spec (A) \ a by gaps of length greater than or equal to d > 0, that is, if 

dist((7,£) =d>0, (1.1) 

then the spectrum of L will also consist of two disjoint components, CO = spec(L) n 
0||V|| (ct) and £2 = spec(L) n 0|| V || (£), provided that 

||V||<d/2. (1.2) 

Under condition (11.2b one may think of the separated spectral components CO and 
£2 of the perturbed operator L as the result of the perturbation of the initial disjoint 
spectral sets a and E, respectively. Clearly, this condition is sharp in the sense that if 
|| V || > d/2, the perturbed operator L may not have separated parts of the spectrum 
at all. 

Assuming ( 11.2b . let P = Ea(<7) and Q — El(co) be the spectral projections of 
the (self-adjoint) operators A and L associated with the unperturbed and perturbed 
isolated spectral sets a and CO, respectively. A still unsolved problem is to answer 
the following fundamental question: 

( i) Is it true that under the single spectral condition ( 11. Il l the perturbation bound 
( 11.2b necessarily implies 

P»-G||<17 (1.3) 

Our guess is that the answer to the question (i) should be positive, but at the moment 
this is only a conjecture. 

Notice that the quantity 9 = arcsin(||P— expresses the maximal angle 
between the ranges Ran(P) and Ran(Q) of the orthogonal projections P and Q 
(see Definition 12. H and Remark l2~2l below) . If 6 < 7C/2, the subspaces Ran(f) and 
Ran((2) are said to be in the acute-angle case. Thus, there is an equivalent geometric 
formulation of the question (i): Does the perturbation bound ( 11.2b together with the 
single spectral condition ( 11.11 ) always imply that the spectral subspaces of A and L 
associated with the respective unperturbed and perturbed spectral sets a and co are 
in the acute-angle case? 

Furthermore, provided it is established that ( 11.3b holds, at least for 

||V||<cd (1.4) 
with some constant c < 1/2, another important question arises: 
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(ii) What function M(x), x G [0, c), is best possible in the bound 
arcsin(||P- Q||) < M l ^-j ? 




(1.5) 



The estimating function M in ( 11.51 ) is required to be universal in the sense that it 
should be the same for all self-adjoint A and V for which the conditions (II . lb and 
Oi l hold. 

Note that if one adds to ( II . II an extra assumption that one of the sets a and E 
lies in a finite or infinite gap of the other set, say, a lies in a gap of E, the answer 
to the question (i) is known to be positive and the optimal function M in the bound 
( 11.51 ) is given by M(x) — j arcsin(2x), x G [0, j). This is the essence of the Davis- 
Kahan sin 2© Theorem ifTTI . For the same particular mutual positions of the spectral 
sets (7 and E the positive answer to the question (i) and complete answers to the 
question (ii) have also been given (under conditions on ||V|| much weaker than ( 11.2b ) 
in the case where the perturbation V is off-diagonal with respect to the partition 
spec (A) = a U E (see the tan 2© Theorem in [11] and the a priori tan© Theorem in 
El ED; cf. the extensions of the tan2© Theorem in lfT2l[T4ll2T1l ). 

As for the general case where no requirements are imposed on the spectral 
sets a and E except for the separation condition ( 11.11 and no assumptions are made 
on the structure of the perturbation V, we are only aware of the partial answers to 
the questions (i) and (ii) found in [16] and [19|. We underline that both [16] and 
lfl9l only treat the case where the unperturbed operator A is bounded. In j 16'] it was 
proven that inequality (11.3b holds true whenever || V | < c KMM d with 



Explicit expressions for the corresponding estimating functions M found in [ 19 1 and 
ifTBI are given below in Remarks 13.61 and [3~8l respectively. The bound of the form 
( 11.51 ) established in 1 19] is stronger than its predecessor in 1 16 1. 

In the present work, the requirement that the operator A should be bounded is 
withdrawn. Introducing the sequence 



we prove that under conditions ( 11.11 ) and dl.2l ) the following estimate of the form 
(TO) holds: 




(1.6) 



In lfl9l the value of c in the bound ( 11.4b ensuring ( 11.31 ) has been raised to 



c ms = \-^2 =0.432332.... 



(1.7) 
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c *-!6 7Z2T7T4 =0.454169.... (1.10) 



with n # (x) = max {n | n € {0} UN,^<i}. The estimate ( 11.8b is sharper than the 
best previously known bound for \\P — Q\\ from [19| (see Remark 1331 for details). 
Furthermore, this estimate implies that ||P— Q|| < 1 and, thus, that the spectral 
subspaces Ran(P) and Ran(2) are in the acute-angle case whenever ||V|| < c*d, 
where the constant 

K 6 - 2% A + 32a: 2 - 32 

(7T 2 +4) 

is larger (and, hence, closer to the desired 1 /2) than the best previous constant ( 11.7b . 

The plan of the paper is as follows. In Section|2]we recall the notion of maxi- 
mal angle between subspaces of a Hilbert space and recollect necessary definitions 
and facts on pairs of subspaces. In Section|3]we extend the best previously known 
subspace perturbation bound (from lfl9l ) to the case where the unperturbed opera- 
tor A is already allowed to be unbounded. Note that the extended bound is later on 
used in the proof of the estimate ( 11.8b . Section [4] is devoted to deriving two new 
estimates for the variation of a spectral subspace that we call the a priori and a pos- 
teriori generic sin20 bounds. These respective estimates involve the maximal angle 
between a reducing subspace of the perturbed operator L and a spectral subspace of 
the unperturbed operator A or vice versa. The principal result of the present work, 
the estimate (11.81 ). is proven in Section |5J (see Theorem 15 ,4b . Under the assumption 
that V ^ 0, the proof is performed by multiply employing the a priori generic sin 29 
bound and using, step by step, the triangle inequality for the maximal angles be- 
tween the corresponding spectral subspaces of the operator A and two consecutive 
intermediate operators Lj =A + tjV, where t j = Xjd/ 1| V|| , j € {0}UN. Finally, in 
Section|6]we apply the bound (11.81 ) to the Schrodinger operator describing a bound- 
edly perturbed A^-dimensional isotropic quantum harmonic oscillator. 

The following notations are used throughout the paper. By a subspace of a 
Hilbert space we always mean a closed linear subset. The identity operator on a 
subspace (or on the whole Hilbert space) is denoted by Im ; if no confusion arises, 
the index *p is often omitted. The Banach space of bounded linear operators from a 
Hilbert space <p to a Hilbert space is denoted by S 0) and by S (<£) if = <p. 
If P is an orthogonal projection in a Hilbert space S) onto the subspace *p, by f x 
we denote the orthogonal projection onto the orthogonal complement *p x := f) G *}5 
of the subspace The notation Er(a) is used for the spectral projection of a self- 
adjoint operator T associated with a Borel set ff C t. By Q r (o), r > 0, we denote 
the closed r-neighbourhood of a in R, i.e. r (a) = {x g E dist(jc, cr) < r}. 



2. Preliminaries 

The main purpose of this section is to recollect relevant facts on a pair of subspaces 
and the maximal angle between them. 

It is well known that if is a Hilbert space then ||P — Q\\ < 1 for any two 
orthogonal projections P and Q in (see, e.g., []] Section 34]). We start with the 
following definition. 
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Definition 2.1. Let *p and be subspaces of the Hilbert space f) and P and Q the 
orthogonal projections in with Ran(P) = *}2 and Ran(Q) = 0. The quantity 

0(«P,0):=arcsm(||P-e||) 
is called the maximal angle between the subspaces *p and 0. 

Remark 2.2. The concept of maximal angle between subspaces is traced back to 
M.G. Krein, M. A. Krasnoselsky, and D.P. Milman 07). Assuming that OP,0) is 
an ordered pair of subspaces with *}3 ^ {0}, they applied the notion of the (relative) 
maximal angle between *}3 and to the number <p in [0, f 1 introduced by 

sin<p(<}3,0) = sup dist(*,0). 

JteCp, ||jr|| = l 

If both q3 ^ {0} and ^ {0} then 

9(q? ) 0) = max{<p(q3,O) )( p(0,q3)} 

(see, e.g., [8, Example 3.5]) and, in general, <p(*£,0) ^ (p{Q,^)- Unlike <p(*}3,0), 
the maximal angle 0(*}3,0) is always symmetric with respect to the interchange of 
the arguments *}3 and 0. Furthermore, 

<p(0,«P) = <p(93,0) = 0(0,qp) whenever | |P- e|| < 1. 

Remark 2.3. The distance function ii(*p,0) = ||P — Q\\ is a natural metric on the 
set S($j) of all subspaces of the Hilbert space fj. Clearly, for the maximal angle 
0(*P,0) to be another metric on &(£)), only the triangle inequality 

0(q*,0) < B(%&) + 0(91,0), for any *p,0,fK 

is needed to be proven. That 0(^3,0) is indeed a metric on §(Sj) has been shown in 
|9l . In Lemma l2. 15l below we will give an alternative proof of this fact. 

Remark 2.4. (<£ x , X ) = 0). This follows from the equalities 

\\P ± -Q ± \\ = \\(I-P)-(I-Q)\\ = \\P-Ql 
where / is the identity operator on fj. 

Definition 2.5. Two subspaces *}3 and of the Hilbert space are said to be in the 
acute-angle case if ^ {0}, ^ {0}, and 0(^,0) < f , that is, if 

||P-G||<1, (2.1) 

where P and Q are the orthogonal projections in f) with Ran(P) = and Ran((2) = 
0. 

Remark 2.6. We recall that the subspaces and are said to be in the acute case 
ifq3H0 x = ( P ± n0 = {O} (cf„ e.g., ifTTI Definition 3.11). The bound \\P-Q\\ < 1 
implies both *pn0 x = {0} and *p x n0 = {0} (see, e.g., fl3] Theorem 2.2]). Hence, 
if the subspaces *p and are in the acute-angle case, they are automatically in the 
acute case. 
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Remark 2.7. It is known (see, e.g., Ifl5l Corollary 3.4]) that inequality d2.ll) holds 
true (and, thus, and are in the acute-angle case) if and only if the subspace 
is the graph of a bounded linear operator X from the subspace *p to its orthogonal 
complement *p x , i.e. 

Q = 5(X):={x®Xx\xe^}. (2.2) 
In such a case the projection Q admits the representation 

/ (/qj +x*xy l [i v +x*x)- l x* \ 
^ \x(i v +x*x)- 1 x(i v +x*x)- 1 x*J y ' 

with respect the orthogonal decomposition f) = < p©*p ± (cf. [15, Remark 3.6]). 
Moreover, under condition d2.ll ) the orthogonal projections P and Q are unitarily 
equivalent. In particular, 

P = U*QU, 
where the unitary operator U is given by 

( (Zp+rx)- 1 ' 2 -x*(i^+xx*)- 1/2 \ nA . 
u ~ \x(/,p+x*x)-V2 (i v±+ xx*)- l/2 )' } 

Remark 2.8. One verifies by inspection that the unitary operator ( 12.41 ) possesses the 
remarkable properties 

U 2 =(Q X -Q)(P J --P) and Ret/>0, (2.5) 

where Re U = \ (U + U*) denotes the real part of U. 

The concept of direct rotation from one subspace in the Hilbert space to an- 
other was suggested by C. Davis in fTO]. The idea of this concept goes back yet to 
B. Sz.-Nagy (see E2 §105]) and T. Kato (see dl Sections 1.4.6 and 1.6.8]). We 
adopt the following definition of direct rotation (see [11 Proposition 3.3]; cf. ETl 
Definition 2.12]). 

Definition 2.9. Let *p and £3 be subspaces of the Hilbert space f). A unitary operator 
S on f) is called the direct rotation from S) to if 

QS = SP, S 2 = {Q ± -Q)(P ± -P), and Re S > 0, (2.6) 

where P and Q are the orthogonal projections in S) such that Ran(P) = *p and 
Ran(g) =0. 

Remark 2.10. If the subspaces *p and are not in the acute case, the direct rotation 
from *p to exists if and only if 

dim(<p n X ) = dim(0 n *P X ) 

(see ifTTl Proposition 3.2]). If it exists, it is not unique. 

Remark 2.11. If the subspaces *p and are in the acute case then there exists a 
unique direct rotation from *p to (see IfTTl Propositions 3.1 and 3.3] or [21 The- 
orem 2.14]). Comparing ( 12.5b with d2.6l ). one concludes that the unitary operator U 
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given by (12.4b represents the unique direct rotation from the subspace *p to the sub- 
space whenever these subspaces are in the acute-angle case. The direct rotation 
U has the extremal property (see iflOl Theorem 7.1]) 

||ff-/j5||= inf \\U-Ifi\\, (2.7) 

f/eU(qj,fl) 

where U(<P, 0) denotes the set of all unitary operators U on $j such that P = U*QU. 
Equality ( 12.7b says that the direct rotation U is norm closest to the identity operator 
among all unitary operators on 5) mapping *}3 onto 0. 

The operator X in the graph representation ( 12.2b is usually called the angular 
operator for the (ordered) pair of the subspaces *p and 0. The usage of this term is 
motivated by the equality (see, e.g., lfl5l ) 

0(«P,0) = arctan^£*X, (2.8) 

where 0(^3,0) denotes the operator angle between the subspaces *}3 and (mea- 
sured relative to the subspace *}3). One verifies by inspection (see, e.g., Ifl5l Corol- 
lary 3.4]) that, in the acute-angle case, 

sin(0(<p,0)) = HP- e|| = J 1 *, 1 '. = sin||0OP,0)||, (2.9) 
which means, in particular, that 

e(«p,jQ) = ||e($p,fi)||. (2.io) 

Furthermore, the lower bound for the spectrum of the real part of the direct rotation 
423} is given by (cf. E] Remark 2.18]) 

min(spec(ReC/)) = 1 =cos(fl(«P,fl)). (2.11) 

To have a more convenient characterization of the distinction between a uni- 
tary operator and the identity one, we recall the notion of the spectral angle. 

Definition 2.12. Let S be a unitary operator. The number 

&(S)= sup | argz| , argz g (-n, n], 

zGspec(S) 

is called the spectral angle of S. 

Remark 2.13. The size of \\S— 1\\ is easily computed in terms of #(S) and vice versa 
(see iPTl Lemma 2.19]). In particular, 

||S-/||=2sin(^V (2.12) 

Furthermore, 

cos#(S) =min(spec(ReS)). (2.13) 

By comparing ( 12.1 lb with ( 12.131 ) and (12.7b with ( 12.121 ) we immediately arrive 
at the following assertion. 

Proposition 2.14. Let *p and be subspaces of a Hilbert space Sj. Assume that ^3 
and are in the acute-angle case. Then: 
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( i) The maximal angle between *}3 and is nothing but the spectral angle of the 
direct rotation U from *P to 0, i.e., 0(^3,0) = &(U). 

( ii) # (S) > 0) for any unitary S on fj mapping onto 0. 

We conclude the present section with a proof of the triangle inequality for the 
maximal angles between subspaces. 

Lemma 2.15. Let <}3, 0, and y\ be three arbitrary subspaces of the Hilbert space 
$j. The following inequality holds: 

0(q8,0)<0(qj,9 £ t) + 0($H,0). (2.14) 

Proof. If 0$J,5K) + 0($H,0) > 7t/2, inequality d214t holds true since 0(^3,0) < 
n/2 by the definition of the maximal angle. 

Suppose that (33,9V) + 0(*K,0) < n/2. In such a case both the pairs 0P,9t) 
and (9t,0) of the argument subspaces are in the acute-angle case. Then there are 
a unique direct rotation U\ from *}3 to 9\ and a unique direct rotation U2 from fH to 
(see Remark l2.1 lb . By [21 , Lemma 2.22], the spectral angle of the product 
S := U2U1 of the unitary operators U\ and U2 satisfies the bound 

<#(£/!) + #(f/ 2 ). (2.15) 

Notice that by Proposition ^. 14l (i) 

= 9{y,m) and $(U 2 ) = 0(9^,0), (2.16) 

because both U\ and U2 are direct rotations. Since Ran([/i |qj) = and Ran(f/2|£j) = 
JH, the unitary operator 5 maps *p onto 0. Hence, #(5) > 0(*P,0) by Proposition 
I2.14l (ii). Combining this with ( 12.15b and (12.16b completes the proof. □ 

3. An extension of the best previously known bound 

In this section we extend the norm estimate on variation of spectral subspaces of a 
bounded self-adjoint operator A under a bounded self-adjoint perturbation V estab- 
lished recently by K. A. Makarov and A. Seelmann in |[T9l to the case where A is 
allowed to be unbounded. 

We begin with recalling the concept of a strong solution to the operator Sylvester 
equation. 

Definition 3.1. Let Ao and Ai be (possibly unbounded) self-adjoint operators on 
the Hilbert spaces Sjq and 5) 1 , respectively, and Y € 23 (S)o , f> 1 ) ■ A bounded operator 
X e 23(ioo,-f)i) is said to be a strong solution to the Sylvester equation 

XA -A l X = Y (3.1) 

if 

Ran(X| Dom(Ao) ) C Dom(Ai) 

and 

XAof - AiXf = Yf for all / e Dom(A ). 

We will use the following well known result on a sharp norm bound for strong 
solutions to operator Sylvester equations (cf. [5, Theorem 4.9 (i)]). 
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Theorem 3.2. Let Aq, Ai, andY be as in Definition \3.1\ If the spectra of Aq and Ai 
are disjoint, i.e. if 

8 := dist(spec(Ao),spec(Ai)) > 0, 

then the Sylvester equation A3. U has a unique strong solution X € 23 (fjo , fj i ). More- 
over, the solution X satisfies the bound 

11*11 < f 0.2) 

Remark3.3. The fact that the constant c in the estimate 5 \\X\\ <c||F|| for the generic 
disposition of the spectra of Aq and Ai is not greater than n/2 goes back to B. Sz.- 
Nagy and A. Strausz (see [23]). The sharpness of the value c — n/2 has been proven 
by R. McEachin ||20l . In its present form the statement is obtained by combining J2 
Theorem 2.7] and (3] Lemma 4.2]. 



The next statement represents nothing but a corollary to Theorem l3.2l 

Proposition 3.4. (cf. [20 1) Let A and B be possibly unbounded self-adjoint opera- 
tors on the Hilbert space fj with the same domain, i.e. Dom(B) = Dom(A). Assume 
that the closure C = B ~A of the symmetric operator B — A is a bounded self-adjoint 
operator on fj. Then for any two Borel sets ffl,ilcl the following inequality holds: 

dist(ffl,fl)||^A(fl>)E«(a)|| <|||C||, (3.3) 

where Ea(co) and Ei(£2) are the spectral projections of A and B associated with the 
sets CO and Q., respectively. 

Proof. Clearly, it suffices to give a proof only for the case where 

O)Cspec(A), £2cspec(B), and dist(a),n) > 0. (3.4) 

Assuming (13.4b . we set P = Ea(@) and Q = Eb(O). The spectral theorem implies 

P/GDom(A)n ( P forany/eDom(A), (3.5) 

ggeDom(fl)n0 forany^eDom(B), (3.6) 

where *p := Ran(P) and := Ran(2) are the spectral subspaces of the operators A 
and B associated with their respective spectral subsets CO and £2. Due to Dom(Z?) = 
Dom(A) the inclusions (13.51 l and ( 13.6b yield 

Ran (Pe| Dom(fi) ) C Dom(A) n tP- (3-7) 

Since P commutes with A, Q commutes with B, P 2 = P, and Q 2 = Q, from (13.51 >- 
dlTI i it follows that 

PQQBQf -PAPPQf = PCQf forany/eDom(B). (3.8) 

Now let Ao, and B& be the parts of the self-adjoint operators A and B associated 
with their spectral subspaces *p = Ran(P) and = Ran(Q)- That is, 

A m =A\<$ with Dom(A ffl ) = < pnDom(A), (3.9) 

Bn =B\q with Dom(B n ) =0nDom(B) (=0nDom(A)). (3.10) 
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Also set X := P\ Q = PQ\q. Taking into account (1377b . d3~9l ), and ( 13.10b we have 
Ran(X| Dom(Bn) ) cDom(A ffl ) 

and then ( 13.8b implies 

XB a f-A a Xf = PCf for any / e Dom(B n ), 

which means that the operator X is a strong solution to the operator Sylvester equa- 
tion 

XB a -A w X = PC\ a . 
To prove (13 ,3b it only remains to notice that 

\\E A (a>)E B (Sl)\\ = p>G|| = pf||, 
spec(A(i)) = o, spec(B a )-ii, \\PC\ Q \\ < \\C\\ 
and then to apply Theorem l3.2l □ 

Theorem 3.5. Given a (possibly unbounded) self-adjoint operator A on the Hilbert 
space $j, assume that a Borel set O C M is an isolated component of the spectrum of 
A, that is, (7 C spec (A) and 

dist((7,E) =d>0, (3.11) 

where E = spec(A) \ O. Assume, in addition, that V is a bounded self-adjoint oper- 
ator on S) such that 

\\V\\<d/2 (3.12) 

and let T(t), t £ [0, 1], be the spectral projection of the self-adjoint operator 

L,=A + tV, Dom(L f ) =Dom(A), (3.13) 

associated with the closed \\V ^-neighborhood 0||v|| (p) of the set O. The projection 
family {r(f)} re [o ij is norm continuous on the interval [0, 1] and 

arcsin(||r(fr) -r(fl)||) < ^log ( ^~2fr]jyjj ) whenever 0<a<A<l. (3.14) 

Proof. Let (O, = spec(L r ) n 0|| V || (cr) and £2, = spec(L ( ) n 0|| V || (E), t £ [0,1]. Since 
A is a self-adjoint operator and L t is given by ( 13.13b . we have 

O^nn, = 0, spec(L r ) = ffliUfl,, 

and, in fact, 

co, c0,||y||((7) and £2 r cO,|| V ||(E), re [0,1]. (3.15) 
Notice that under condition ( 13.121 ) from ( 13.111 ) and ( 13.151 ) it follows that, for s,t £ 

[0,1], 

dist(a*,n,)></-f||V||-.s||V|| =rf-||V||(r + j). (3.16) 

In particular, 

d-\\V\\(t + s)>d-2\\V\\>0 whenever s,t £ [0,1]. (3.17) 
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Obviously, T(t) = E Ll (co,) and T(t ) x = E L , (Q,), f e [0, 1]. Thus, for any s,t € 
[0, 1] Proposition l3.4l implies 



l|rwr(0 x ll < - mit ~\ and ||r(,) x r (f) || < n - (3.18) 

11 w w 11 " 2 dist(co. s ,£2 ( ) " 2 dist(a,co,) 

Since ||r(f) ~r(j)|| = max{ \\T(t )r(i) x || , ||r(f ) x r(s) || }, from (l3~T6l l and (EH one 
concludes that 



™- r ^''4 ^Fff^ ^ any [0,1]. (3.19) 

In view of ( 13.17b . the operator norm continuity of the projection path {r(f)} fe [o,i] 
on the interval [0, 1] follows immediately from estimate (13.19b . 

Now suppose that s < t (as before, s, t g [0, 1]) and observe that 

t — s [' dr 

< / , (3-20) 



d-\\V\\(t + s) J s d-2\\V\\x 

Indeed, the difference between the right-hand side and left-hand side parts of (13.20b 
may be written as 

| TC (/(t)+/(2t c - t) -2f(Xc))dz, (3.21) 
where T c = (s + t)/2 is the center of the interval [s,t] and 

One verifies by inspection that the expression /(t) +/(2t c — t) — 2f(x c ) under 
the integration sign in ( 13.21b is strictly positive for T £ [s, T c ) and zero for T = T c . 
Therefore, the integral in (13.21b is positive and hence inequality ( 13.20b holds true. 

Assume that < a < b < 1. For a sequence of points fo,?i, • • • ,t n € [a,b], n € N, 
such that 

a = t <h <...<t n = b (3.22) 
by d3~T9b and d3~20l ) one obtains 

I l|r(o +1 -r(»y)|| E 

ttII VII 'i^ 1 /-o+i rfi 



2 £ 4 d-2||V||T 
7r||V|| [ b dx 



2 y fl d-2\\v\\x (123) 

Evaluating the last integral in ( 13.23b and taking supremum over all choices of n e N 
and fo,fi , . . . , t„ £ [a,b] satisfying (13.22b results in the bound 

. ( r)<?io g ^ 2fl " y " 



4 6 v^-2^iivii;' 
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where 



n-l 



^(r):=su P £ \\r(t J+1 )-r(tj) 



n e N, a = t < t\ < . . . < t n = b > (3.24) 



is the length of the (continuous) projection path T(t ), f G [a,b]. Applying lfl9l Corol- 
lary 4.2], which establishes that arcsin(||r(fe) — r(a)||) < completes the proof. 

□ 

Remark 3.6. Let 21 = Ran(E A (cr)) and £ = Ran(E £ ((o)) where L := A + V with 
Dom(L) = Dom(A) and co = spec(L) n 0||yii(a). By setting a = and b = 1 in 
( I3.14l i. one obtains 

0(21,£)<M MS (W) ; ( 3. 25) 

where (21, £) is the maximal angle between the spectral subspaces 21 and £ and 

M Ms(*) : =f lo s(_L;) > , 6 [0,i). (3.26) 

For bounded A, the estimate (I3.25l l has been established in [19, Theorem 6.1]. 

Note that ( 13.25b implies that 0(2l,£) < f and, thus, that the subspaces 2t and 
£ are in the acute-angle case if 

\\V\\<c MS d, (3.27) 

where c MS is given by (11.7b . 

Remark 3.7. For future references we remark that, due to (13.25b . 

0(2t,£)<! whenever ||V|| < c %/A d, (3.28) 

where 

c ,=]-- — = 0.3 16060.... (3.29) 
2 2e 

Remark 3.8. The bound ( 13.25b is stronger than the earlier estimate ifTBI 

(<f 

where the value of c,^,^ is given by ( 11.6b and 



0(2t,£)<M KMM M-l <- if||F||< CKMM rf, (3.30) 



/ fly \ 

m kmmW— arcsinf 2(1 _ x) J, 0<x<c mM . (3.31) 
The estimate (13.30b was established in the proof of Lemma 2.2 in IfTBI . 



4. A priori and a posteriori generic sin 20 estimates 

We begin this section with the proof of an estimate for sin(20(2l,£)), where 21 is 
a reducing subspace of the self-adjoint operator A and £ is the spectral subspace 
of the boundedly perturbed self-adjoint operator L=A + V. In general, 21 does not 
need to be a spectral subspace of A. 
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Theorem 4.1. Let A be a (possibly unbounded) self-adjoint operator on the Hilbert 
space 9). Suppose thatV is a bounded self-adjoint operator on Sj andL = A +V with 
Dom(L) = Dom(A). Assume that 21 is a reducing subspace of A and £ is a spectral 
subspace of L associated with a Borel subset CO of its spectrum. If the subspaces 21 
and £ are in the acute-angle case then 

dist((»,fl)sin(20) <fc||V||, (4.1) 

where D. — spec(L) \ co is the remainder of the spectrum of L and := 0(2l,£) 
denotes the maximal angle between 21 and £. 

Proof. For dist(o,n) = the bound ( 14. U is trivial. Throughout the proof below 
we will assume that dist(a),£2) ^ 0. 

Since 21 is a reducing subspace of the self-adjoint operator A, its orthogonal 
complement 21 x = 21 is also a reducing subspace of A. Furthermore, Dom(A) = 
Dom(Ao) ©Dom(Ai), where Dom(Ao) and Dom(Ai) are domains of the parts Aq = 
A|«a and Ai = A| a ± of A in its reducing subspaces 21 and 2l x , respectively (see, e.g., 
Section 3.6]). 

Let P be the orthogonal projection in Sj onto the subspace 21. Since V G 23 (io), 
the (self-adjoint) operator L admits the following block representation with respect 
to the orthogonal decomposition fj = 21© 2t x : 

L= (b* Di )' Dom ( L )= Dom ( z) o)©Dom(D 1 ) (=Dom(A)), (4.2) 

where B = PV| a ±, D = A +PV | a , with Dom(D ) = Dom(A ) and D l = A x + 
P^Vl^x, withDom(Di) =Dom(A!). 

That the subspaces 21 and £ are in the acute-angle case implies that there is 
a bounded operator X from 21 to 2l x such that £ is the graph of X, that is, £ = 
SPO (see Remark |2~7] |. It is well known (see, e.g., J2] Lemma 5.3]) that the graph 
subspace SPO is a reducing subspace for the block operator matrix ( 14. 2t if and only 
if the angular operator X is a strong solution to the operator Riccati equation 

XDq-DxX+XBX =B*. (4.3) 

The notion of the strong solution to (14.3b means that (see [2 3 |; cf. Definition [3Tj 

Ran(X| Dom(Do) ) cDom(J)i) (4.4) 

and 

XD f-D l Xf + XBXf = B* f for all /eDom(D ). (4.5) 

It is straightforward to verify that if X is a strong solution to ( 14.5b then 

XZof-ZiXf = B*(I + X*X)f forall /eDom(Do), (4.6) 

where Z = D Q +BX with Dom(Z ) = Dom(D ) = Dom(A ) and Z x = D x -B*X* 
with Dom(Z!) = Dom(D 1 ) = Dom(Ai). 

Our next step is in transforming (14.6b into 

x(i +x*xy 1/2 A (i +x*x) l i 2 f - (i +xx*y 1/2 A l (i +xx*) l/2 xf 

= B*(I + X*X)f forall/eDom(D ), (4.7) 
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where Ao and Ai are given by 

Ao = (I + X*X) l / 2 Z Q (I + X*X)- 1 / 2 , 

Dom(Ao) =Ran((/ + X*X) I /2| Dom(Do) ), (4 ' 8) 

and 

Ai = (/ + XX*) l / 2 Zi(I + XX*y l l 2 , 

Dom(A 1 ) =R a n((I + XX*y/ 2 \ Dom(Di) ). (4 ' 9) 

Note that the (self-adjoint) operator L is unitary equivalent to the block diagonal 
operator matrix A = diag(Ao,Ai), Dom(A) = Dom(Ao) ©Dom(Ai) (see, e.g., J2] 
Theorem 5.5]) and, thus, both Ao and Ai are self-adjoint operators. Furthermore, 
from ( 14. 8t and ( 14.9t it follows that Zq and Z\ are similar to Ao and Ai, respectively. 
Combining [2, Theorem 5.5] with J3] Corollary 2.9 (ii)] then yields 

spec(Ao) = spec(L|£) = co and spec(Ai) = spec(L| £ ±) = Q.. (4.10) 

Applying (I + XX*) l l 2 from the left to both sides of ( 14.6b and choosing / = 
(/ + X*X)~ x l 2 g with g € Dom(Ao), we arrive at the Sylvester equation 

KA g-A 1 Kg = Yg forallgeDom(Ao), (4.11) 

where 

K= (I + XX*y/ 2 X(I + X*X)- 1/2 , (4.12) 

Y =(I + XX*) 1/2 B*{I + X*X) 1/2 . (4.13) 

we have Ran( (I+X*X)~ l / 2 \ !A) )= Dom(D ). Furthermore, 



Byc^.w.a.wvo.n^ iDom(Ao) 

Ran(X| Dom (D )) cDom(Di) 

by d4~4] i, and thus, by ( FT91 . 

Ran (^lDo m( A ))cDom(A 1 ). (4.14) 

Hence K is a strong solution to the Sylvester equation ( 14.1 11 1. 

It is easy to verify (see 161 Lemma 2.5]) that {I+XX*) l l 2 X =X(I +X*X) l l 2 . 
Thus, (14.12b simplifies to nothing but the identity X = K, which by (14. 1U and ( 14. 14t 
means that X is a strong solution to the Sylvester equation 

XKq-KiX = Y. (4.15) 

Observe that ||F|| < ||B|| (1 + ||X|| 2 ). Taking into account ( BTTOt , Theoreml3~2lvields 

11X11 < % \\B\\ 



1 + ||X|| 2 - 2 dist(fi>,Q)' 
Now the claim follows from the fact that 2||X||/(1 + ||X|| 2 ) = sin(20) by dPt , □ 

Remark 4.2. Clearly, if 7r||V|| > dist(a),i2), the estimate ( 14. \ \ is of no interest. Sup- 
pose that dist(a),i2) = 8 > and 7l||V|| < 8. In such a case ( 14. Il l does allow to 
obtain a bound for the maximal angle 9 but only provided the location of 6 relative 
to j is known. In particular, if it is known that 9 < j then ( 14.11 ) implies the upper 
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On the contrary, if it is known that 9 > ? then ( 14. It 



yields the lower bound > § — j arcsin ^^5-^ 

Corollary 4.3. Lef A foe a (possibly unbounded) self-adjoint operator on the Hilbert 
space Sj. Assume that 21 is the spectral subspace of A associated with a Borel subset 
O of its spectrum. Suppose that V is a bounded self-adjoint operator on and 
L = A + V with Dom(L) = Dom(A). Furthermore, assume that £ is a reducing 
subspace ofL. If the subspaces 21 and £ are in the acute-angle case then 

dist(ff,£)sin(20)<jf||V||, (4.16) 

where E = spec(A) \ <J is the remainder of the spectrum of A and 9 denotes the 
maximal angle between 21 and £. 

Proof. Consider A as the perturbation of the operator L, namely view A as A = L + W 
with W = —V, and then the assertion follows from Theorem l4.ll □ 



Remark 4.4. Suppose that dist(a,E) = d > and ||V|| < c.^d, where c , 4 is 
given by d3~29l ; observe that c n/4 < \ = 0.318309... and, thus, ^ < 1. Let 
CO = spec(L) n 1| v || (c) and £ = Ran(Ei(o))). By Remark l3~7l under the condition 
\\V\\ < c, 4 d we have 9 ('01,2) < f and, hence, (I4.16l l yields the bound (cf. Remark 

1 ^7t\\v\\\ 

0(2l,£) < - arcsin f J . (4.17) 

Notice that 

1 0.288400... and c,, > -5— r- (4.18) 



7T 2 + 4 ' "' "I 4 K 2 + A 

In Theorem 15 .41 below we will prove that for -^d < ||V|| < \d there is a bound 
on 0(2l,£) tighter than estimate (14.17b . 

The estimates (14. Il l and ( I4.16l l we obtained in Theorem |4T| and Corollary 14. 31 
will be called the generic a posteriori and a priori sin 29 bounds, respectively. These 
estimates resemble the corresponding bounds from the celebrated sin 2© theorems 
proven by C. Davis and W. M. Kahan in ifTTl Section 7] for particular dispositions 
of the sets CO and Q. or a and E. Recall that, when proving the sin 2© theorems, it is 
assumed in IfTTl that the convex hull of one of the sets CO and D. (resp., the convex 
hull of one of the sets a and E) does not intersect the other set. An immediately 
visible difference is that the constant % shows up on the right-hand side parts of 
the generic sin 29 estimates (14. jj and (14.161 1 instead of the constant 2 in the Davis- 
Kahan sin2@ theorems. 

5. A new rotation bound by multiply employing sin 29 estimate 

Let {x,i}n = o be a number sequence with 

4 7t 2 -4 

^o=0, x n = + A ^n-i, n=l,2,.... (5.1) 

7T z +4 7T z +4 
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One easily verifies that the general term of this sequence reads 

l-q" 



11 = 0,1,2,..., (5.2) 



9 = ^-nr< 1 - ( 5 - 3 ) 



2 ' 

with 

n 2 -4 
n 2 + 4- 

Obviously, the sequence ( 15.2b is strictly monotonously increasing and it is bounded 
from above by 1/2, 

x < hi < x 2 < . . . < x„ < . . . < 1/2. (5.4) 

Moreover, 

limx„ = l/2. (5.5) 

n—^oa 

Thus, this sequence produces a countable partition of the interval [0, j) 9 

[0,i) = Q[x B ,x n+ i). (5.6) 

n=0 

Taking into account ( 15.6b . with the sequence {x„}~ =0 we associate a func- 
tion M*(jc), M* : [0, |) — > M, that is defined separately on each elementary interval 

[x B ,x„+i) by 

, . I n / 4tt \ 1 (nix — x„)\ 

^W kK,,)^™ 1 o„ h b = 0,1,2,.... 



t^^n+i) 2 V^+V 2 V l-2x„ 

(5.7) 

We note that the function M*(x) may be equivalently written in the form jl.9\ . 

Proposition 5.1. The function M*{x) is continuous and continuously differentiable 
on the interval [0, 5). Furthermore, this function is strictly monotonously increasing 
on [0, |) ant/ limM*(x) = +°°. 

Proof. Clearly, one needs to check continuity and continuous differentiability of 

M+{x) only at the points x„, n = 1,2, Given n € N, by (15.7b for x 6 [x„_i,x„) 

we have 



M,(x) = — arcsm j + -arcsm ^ j-^— ) , * e 



(5.8) 

Note that (|5~TT > yields 

4 a: 2 -4 

x„ - x„_i = + - 

= ^5— -(l-2x„_i), forallneN. (5.9) 

By ( 15. 9t one observes that 

■r-Xn-i ^ x„ -x„_i _ 4 ^ ^ 

1 — 2x„_ 1 1 — 2x„_ 1 K 2 + 4 
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and then from ( 15.8b it follows that 

n ( An \ 
maMJx) = -arcsin . (5.11) 

Meanwhile, by its definition ( 15.7b on the interval x G [xr„, >f„+i), the function M+(x) 
is right-continuous at x = x n and 

n / \% \ . . , . 

-arcsin 7 ■ - = hm M*(jc) = M*(x„). 
2 \ K z + 4 / x4j«j, 

Hence, ( 15.1 11 1 yields continuity ofM*(x) at * = k„. 

Using equality d5.9l l one more time, one verifies by inspection that for any 
n EN the left and right limiting values 

lim<(x) = - 1 === ^ and \imM',{x) = \ % 



1 



n\>c n - >c n _ 

(1 



of the derivative M+{x) as x — > x„ are equal to each other. Hence, for any n G N the 
derivative M'^x) is continuous at x = x n and then it is continuous on [0, A). 

Obviously, M'+(x) > for all x £ (0, |), which means that the function M+(x) 
is strictly monotonously increasing on [0, |) . Then, given any n € N, from ( 15.7b it 

follows that M*(x) > nCo with Co = \ arcsin ( J*^ 4 J > whenever x > x n . Taking 

into account ( 15.5b , this implies thatM*(x) —> °° as x —> A, completing the proof. □ 

Remark 5.2. Since the function M+(x) is continuous and strictly monotonous on 
[0, A) and M*(0) = 0, limM*(x) = +°°, there is a unique number G [0, A) such 

that 

M*(c*) = |. (5.12) 
An explicit numerical evaluation of 

1 / 4n \ 

C = - arcsin ) (5-13) 

shows that C = § x 0.360907.... Thus, 2C < f while 3C > f . Hence, c* G 
[>T2j >q) and by ( 15.7b . with n = 2, equation ( 15.12b turns into 

anaBn l^T4 j + 2 " rc " m J = 2 ' (5 ' 14) 

with ^ being equal (see (15.2b ) to 

x 2 = . - 87t = 0.410451 ... f>-). (5.15) 
(7T 2 + 4) 2 V ^7 

One verifies by inspection that the solution to equation ( 15.141 ). and hence to equa- 
tion ( 15.121 ). reads 

7T 6 -27r 4 + 327r 2 -32 
c* = 16 , 7 -r-, =0.454169.... (5.16) 

[n 2 + A) 4 



18 



S. Albeverio and A. K. Motovilov 



Proposition 5.3. (Optimality of the function M*.) Let {jU«}"=o c [P> l) ^ e an arbi- 
trary monotonously increasing number sequence such that 

Mo = and 0< ^ ~ < t f or n >\, (5.17) 

1 - 2^„_ 1 

Assume that sup n \l n := /x sup ana! introduce the function F : [0,/X sup ) — > K fey 

F W|[o. M i) = 2 arcsin ( 7r ' x: )' (5 - 18) 

(5.19) 

The function M*(x) is optimal in the sense that if {/^n}" = o afoes nof coincide with 
the sequence (15. Il l, fnen fnere a/ways erafs an open interval 5" C (0,/i sup ) SMcn fnaf 
F(x) > M*(jt) /or a/Z x £ J. 

Proof. First, we remark that if one chooses jj.„ = tt n , n — 0, 1,2 . . . , where {>%}" = q 
is the sequence ( 15. Il l, then F coincides with M*. If {£in}™=o 7^ { x n}~=o then there 
ismSN such that /i„ = x„ for all n < m and fi m 7^ Since, Mm-i = x„,_i, there 
are two options: either 

x m -i< \l m < min{ k,„ , fi m+ 1 } (5 .20) 

or 

x m < jJ, m (and jx m < jl m+ \). (5.21) 

By (E2J and (EH), (l5~T9l in the case (l5T20t equality M*(x) = F(x) holds for all 
x e [0,li m ]. Forx e (n m ,nan{x m ,fi m+1 }) we have 

M*(x) = (m - 1 Co + - arcsin ) - U 
F(x) = (m— 1)Cq + - arcsin I + - arcsin 1 



2 V 1 — 2>£ m _i 12 V 1 — 2jti„ 

where Co is given by ( 15.13b . Having explicitly computed the derivatives of M*(x) 
and F(x) for x G (/ii ,min{^ m ,/x,„+i), one obtains that F'(x) > M^(x) whenever 
\l m < x < % m , where 

§ m = minj^+j , [4 + (7r 2 - 4) (x m _ 1 + Mm)] } . 

Notice that § m > \i m since Mm+i > \l m and 

1 , 7T 2 +4/ 4 7T 2 -4 

^[4 + (* -4)(* m -t + Mm)] - Mm = -^s- + ^^-1 -M« 

— 2^-2 
>0 

(see ( 15. U and ( 15.201 )). Observing that M*(n m ) =F(p. m ), one concludes that F(x) > 
M- k (x) at least for all x from the open interval 3^ = (jj, m ,^ m ). 
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If (15.211 1 holds then M* (x) = F (x) for x £ [0, x m ] . At the same time, for 

K m <x< min{ x m+u ji m } (5 .22) 

we have 

M*{x) = (m - 1)C + - arcsin -A — i + - arcsin -f— - '- 

L \ 1— 2x m _i / 2 \ 1 — 2x m 

(5.23) 

F(x) = (m- l)C + ~ arcsin . (5.24) 

One verifies by inspection that under condition ( 15.221 ) the requirement F'(x) > 
(x) is equivalent to 

1 i 
x > ^ [4 + - 4) (x„,_ i + x„,)] 

_7T 2 +4/ 4 7T 2 -4 7T 2 -4 

7t 2 +4 ( 7t 2 -4 



2n 2 \ m k 

= x,„, 

by taking into account relations ( 15. U at the second step. That is, ( 15.221 ) implies 
F'(x) >M[(x). FromM*(x m ) = F(x m ) it then follows that F(jc) >M*(x) at least for 
all x from the open interval "J = (x„,,min{x m+ i, /!„,}). The proof is complete. □ 

Finally, we turn to the main result of this work. 

Theorem 5.4. Given a (possibly unbounded) self-adjoint operator A on the Hilbert 
space Sj, assume that a Borel set O C R is an isolated component of the spectrum of 
A, i.e. a C spec (A) and 

dist(<7,£) = d > 0, 

where E = spec (A) \ <J. Let V be a bounded self-adjoint operator on Sj such that 

\\V\\<d/2 (5.25) 

and let L = A + V with Dom(L) = Dom(A). Then the maximal angle 9 (21, £) be- 
tween the spectral subspaces 21 = Ran(E,4(<7)) and £ = Ran(Ei(o))) of A and L 
associated with their respective spectral subsets a and CO = spec(L) CI CD || V|| ( c7 ) sat ~ 
isfies the bound 

d(%,£)<M*n^J, (5.26) 

where the estimating function M^{x), x € [0, 5), is defined by (15.7I ). In particular, if 

\\V\\<c,d, (5.27) 

where c* is given by ( 15.161 ), then the subspaces 21 and £ are in the acute-angle case, 
i.e. 0(2l,£) < f. 
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Proof. Throughout the proof we assume that ||V|| ^ and set 

IIVII 

x=±-±. (5.28) 
a 

The assumption ( 15.25b implies x < j . Hence, there is a number n e NU{0} such that 
x £ [x„,x n +i) with x„ and x n+ \ the consecutive elements of the sequence ( 15.1b . 

For n = the bound ( 15.261 ) holds by Remark 14.41 since in this case < 
^4 < c , 4 and M*(x) = ± arcsin(7Tx) (see (14.18b . (l5~TT l. and (E2J). 

In the case where n > 1 we introduce the operators 

Vj = Kj '\^\\ V and + Dom ( L >) = Dom(A), 7 = 0,1,..., re, 

where xy are elements of the sequence (15.1b . Since Xj < x n for j < n and x n < x, 
from ( 15.28b it follows that 

||V>|| = -^||7||<||7||<-, ;=0,l,...,re. 

Therefore, the spectrum of the (self-adjoint) operator Lj consists of the two disjoint 
components 

toy = spec(L ; -) n 0||y.|| (a) and Q.j = spec(Ly) n 0||y|| (E), = 0, 1, ...,re. 
Moreover, 

5;:=dist(ffly,ay) >d-2||7y|| = d(l-2x/), y = 0,l,...,/i. (5.29) 

By we will denote the spectral subspace of Lj associated with its spectral com- 
ponent coj, i.e. = Ran(E Lj (t»y)) . Notice that L = A, coq = a, and £ = 2L 

For < j < re — 1 the operator L; + \ may be viewed as a perturbation of the 
operator Lj, namely 

L j+l= L j+Wj, 7=0,l,...,re-l, 

where 

W;:=V; +1 -V; = (x; +1 -^)^V. 

Similarly, we write 

l = l„+w, with w:=y-y„ = (x-x„)^v. 

One easily verifies that 

C0j + i =spec(L /+1 )nO|| W .|| (a) ; -), j = 0,l,...,«-l, (5.30) 
0) = spec(L)nO|| W ||(a)„). (5.31) 
By taking into account first (15.29b and then ( 15.9b , one observes that 

M = fe^<^^= * j = ^_ n _ h (5 . 32) 

Oj Oj l—2Xj 7T z +4 

and 

\\W\\ _ (x-x„)d x-x„ f x n+ i - x n _ 4 \ 

Sn ~ 8 n ~ 1-2X„ \ 1-2X„ ~ 7t 2 +4 r ( ' 
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0.1 0.2 0.3 0.4 x 



FIGURE 1. Graphs of the functions ^M KMM (x), |M MS (x), and 
^iM+ix) while their values do not exceed 1. The upper curve de- 
picts the graph of -|Mkmm (x) for x < c KMM , the intermediate curve 
represents the graph of |MmsM forx < c MS , and the lower curve 
is the graph of ^M*(x) forx < c*. 



Recall that ^ < c k/4 (see (EH and KWi ). Thus, by ( EHOl l-l EHl RemarkEy] 
applies to any of the pairs (£ ; -,£y + i), j = 0, 1,. ..,n— 1, and (£„,£), which means 
that 

1 . /7r||W,||\ ,1 . / All 



+ 1 



< 



9 (£„,£) < -arcsin 



7r||W/ 



1 

< - arcsin 
~ 2 



< - arcsin 
~ 2 



2 +4 
1 - 2x„ 



7 = 0,1,. 



,n—l, 
(5.34) 

(5.35) 



Meanwhile, by using the triangle inequality for maximal angles between subspaces 
(see Lemma l2.15b one obtains, step by step, 

0(£o,£)<0(£o,£i) + 0(£i,£) (if«>l) 

< 0(£ o ,£i) + 0(£i,£ 2 ) + 0(£2,£) (if«>2) 

< ... 



n-\ 



(ifn > 1) 



(5.36) 



< £fl(£y,£y+i) + ©(£«,£) 

Combining ( 15.36b with ( 15.34b and (15.35b results just in the bound (15.26b . taking into 
account equality £o = 21 and the definitions (15.28b of x and (15.7b of M*(x). Finally, 
by Proposition 15 . 1 1 and Remark IB~2l one concludes that under condition ( 15.27b the 
subspaces 21 and £ are in the acute-angle case. This completes the proof. □ 
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Remark 5.5. Recall that the previously known estimating functions for 0(21, £) are 
those of references lfl6l and fl9l . namely the functions M KMM (^f-^j (see Remark 

13 .8b andM MS ^^j^ (see Remark [3~6l l. One verifies by inspection that the derivatives 

M LmM' M LsM' and M *M of the estimatin g functions M mM (x), M MS (x), and 
M+(x) possess the following properties 

K S ( x )< M Lm( x ) forall^e(0,c KMM ), (5.37) 

M *M < Ks (*) for a11 x e [°. j) \ W"=o> (5.38) 
M:(x„)=M; s (x„), n = 0,l,2,..., (5.39) 

where c KMM is given by ( 11.6b and {x n }~ = o is the sequence (15.1b . Since 

M KMM (0)=M MS (0)=M*(0)=0, 
from (EEJ-fEU it follows that 

M*(x) < M MS (x) < M KMM {x) for all x £ (O.c^J 

and 

M* (*) < M MS (x) for all * <E [c KMM , i ) . 

Thus, the bound (15.26b is stronger than both the previously known bounds |[T6l [T9l 
for 0(2l,£), in particular it is stronger than the best of them, the bound ( 13.251 ). 
established in lfl9l . 

For convenience of the reader, the graphs of the estimating functions M*(jc), 
M MS (x), and M KMM (x), all the three divided by n/2, are plotted in Fig. Q] Plotting of 
the function §M KMM (x) is naturally restricted to its domain [0,c KMM ]. The functions 
M*(x) andM MS (x) are plotted respectively for x £ [0,c*] andxG [0,c MS ] where c* is 
given by (15.16b and c MS by (11.7b . 

We conclude this section with an a posteriori result that is an immediate corol- 
lary to Theorem l5.4l 

Theorem 5.6. Assume that A and V are self-adjoint operators on the Hilbert space 
fj. Let V £ andL—A + V withDom(L) = Dom (A). Assume, in addition, that 

CO is an isolated component of the spectrum of L, i.e. dist(<B,i2) — 8 > 0, where 
D. = spec(L) \ CO, and suppose that \\V\\ < 8/2. Then the maximal angle 0(2l,£) 
between the spectral subspaces 21 = Ran(EA(cT)) and £ = Ran(Ei(o)) of A and 
L associated with their respective spectral components O = spec(A) n 0||v|| (®) an d 
CO satisfies the bound 

0(a,£)<M*fEi^ (5.40) 

with M* given by ( 15.7b . In particular, if\\V\\ < c+8 where c+ is given by (15.16b , the 
subspaces 21 and £ are in the acute-angle case. 

Proof. Do exactly the same step as we did in the proof of Corollary 14.31 Represent 
A as A = L + W with W — — V. Then the assertion follows from Theorem l5.4l □ 
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Remark 5.7. As usually, let E = spec(A) \ a and d = dist(a,E). Suppose that both 
the 'a priori' and 'a posteriori' distances d and 8 are known. Then, depending on 
which of the distances d and 8 is larger, one may choose among the bounds ( 15.26b 
and ( 15.401 ) the stronger one: 

9(%£)<mJ y 

\max{a, 8} ) 



6. Quantum harmonic oscillator under a bounded perturbation 

In this section we apply the results of the previous section to the A^-dimensional 
isotropic quantum harmonic oscillator under a bounded self-adjoint perturbation. 

Let = L2(S. N ) for some N E N. Under the assumption that the units are 
chosen such that the reduced Planck constant, mass of the particle, and the angular 
frequency are all equal to one, the Hamiltonian of the isotropic quantum harmonic 
oscillator is given by 

(A/)(x) = -lA/(x) + i|x| 2 /M, 

(6.1) 



Dom(A) = <f€W : 



dx\x\ 4 \f(x)\ 2 <. 



where A is the Laplacian and W^(M. N ) denotes the Sobolev space ofL 2 (R A ')-functions 
that have their second partial derivatives in L2(M. N ). 

It is well known that the Hamiltonian A is a self-adjoint operator in L2(M. N ). 
Its spectrum consists of eigenvalues of the form 

k n = n+N/2, ra = 0,l,2,..., (6.2) 

whose multiplicities m n are given by the binomial coefficients (see, e.g., |[T8| and 
references therein) 



N + n- 1 
n 



0.1.2 (6.3) 



For n even, the corresponding eigenfunctions f(x) are symmetric with respect to 
space reflection* M> —x (i.e. f(—x) — f(x)). For n odd, the eigenfunctions are anti- 
symmetric (i.e. f(—x) = —f(x)). Thus, if one partitions the spectrum spec(A) = 
cr U X into the two parts 

a = {n+N/2 | n = 0,2,4,...} and L = {n +N/2 \ n= 1,3,5...}, 

then the complementary subspaces 

2l=L 2 , even (E iV ), a J -=L 2 ^ jdd (R JV ) (6.4) 

of symmetric and anti-symmetric functions are the spectral subspaces of A corre- 
sponding to the spectral components a and E, respectively. Clearly, 

<i = dist((7,E) = 1. 

Let V be an arbitrary bounded self-adjoint operator on L2(M. N ) such that 

< 1/2. (6.5) 
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The perturbed oscillator Hamiltonian L = A + V , Dom(L) = Dom(A), is self-adjoint 
and its spectrum remains discrete. Moreover, the closed 1 1 V | -neighborhood | ^ (A„) 
of the eigenvalue (16.2b of A contains exactly m„ eigenvalues X' k , k = 1,2, . . . , m n , 
of L, counted with multiplicities, where m„ is given by (16.3b (see, e.g., |[T3l Section 



where is given by (15.161 1. Let £ be the spectral subspace of the perturbed Hamil- 
tonian L associated with the spectral subset co = spec(L) n 0||y|| (a). Theorem 15. 41 
ensures that under condition ( 16.6l l the unperturbed and perturbed spectral subspaces 
21 and £ are in the acute-angle case. Moreover, the maximal angle 0(2t,£) between 
these subspaces satisfies the bound 



where stands for the function given by (15.7b . 

Obviously, under condition (16.6b the orthogonal complement £ x is the spectral 
subspace of L associated with the spectral set £2 = spec(L) n 0||y|| (£). By Remark 
12.41 0(2l x ,£ x ) = 0(2l,£). Hence the maximal angle between the subspaces 21 x 
and £ x satisfies the same bound ( RT7l i. i.e. d^,^) < M*(||Vj|). Surely, this can 
also be seen directly from Theorem 15. 41 
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